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PART I— THEORY 

( Abstract - -Long Island simulates in a general way an aquifer in the form of an infinite strip 
confined between parallel boundaries at constant head (sea-level), over which recharge from pre- 
cipitation is assumedly uniform. The non-steady flow of water in this idealized system is analyzed, 
assuming provisionally that the effective thickness of saturated beds below sea-level is great com- 
pared to the maximum height of the water-table above sea-level. The rate of accretion to the 
water-table is assumed to vary discontinuously, supposedly being constant for each of the succes- 
sive periods (yearly or monthly) and proportional to the average rate of precipitation during that 
period. The decay of the water-table profile, beginning with any one of the succession of super- 
posed non-steady states, is shown to follow in general a relation composed of terms varying with 
time as exp(-t/tQ), in which tg is a function of the effective porosity, the thickness and the trans- 
mission-constant of the aquifer. This exponential curve may be approximated by a parabola which 
is used to determine values of “effective average rate of precipitation” from published records of 
annual or monthly precipitation. By the “effective average rate of precipitation” at any time is 
meant that rate of precipitation which, had it been maintained uninterruptedly throughout the past, 
would have produced the same water-table profile as actually existed at that particular time. It 
is demonstrated that the effective average rate of precipitation may be determined also simply by 
cumulating departures from progressive averages of precipitation, multiplying the values thus de- 
termined by a known rational coefficient, and adding the appropriate initial value of effective aver- 
age precipitation.) 

564 
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In attempts to correlate precipitation-data with records of ground-water level or stream-flow 
it has long been the custom to cumulate departures of the precipitation from the “normal’' or 
mean. So-called “residual mass-curves” have been drawn showing the excess or deficiency of 
rainfall over the mean for a fixed period. Early work along this line was done in connection with 
the Wiener Neustadt Project for the well-water supply of Vienna [see 1 of “References” at end of 
paper]. 

More recently progressive averages have also been used. For example, a publication on rain- 
fall and runoff [2] gives graphs showing the “annual and ten-year progressive average precipitation” 
for selected drainage-basins in the United States. In most cases on visual comparison the ten-year 
progressive average precipitation shows general conformity with the ten-year progressive average 
runoff. 

It is the intent of this paper to justify use of the cumulative departure from progressive aver - 
ages for correlating precipitation with ground- water levels and with “ground- water flow” in 
streams, thus combining in a sense the idea of cumulating departures with the idea of using pro- 
gressive averages. To the writer’s knowledge, the first attempts along this line in hydrology were 
made by LEGGETTE [3]. In studying fluctuations of the ground-water level on Long Island, New 
York, Leggette cumulated departures of precipitation from ten-year progressive monthly averages. 
That is, he took the difference between the precipitation during a given month (average for two sta- 
tions- -New York City and Setauket, Long Island) and the average precipitation for the ten preceding 
months of the same name, cumulated such differences, and compared the resulting graph with an 
average water-level graph [4]. With the scales of plotting arbitrarily fixed, general overall con- 
formity between the two graphs is suggested for the last nine years (1932-40). Leggette states 
that “a comparison between cumulative departures based on three-year progressive averages and 
cumulative departures based on ten-year progressive averages indicates that the ten-year pro- 
gressive averages are somewhat better for showing the relation between precipitation and ground- 
water levels.” As discussed in Part II of this paper [see abstract at end of Part I], it has since 
been foxmd that cumulative departures based on 25-year progressive averages correlate best with 
ground-water levels on Long Island. 

This paper is an extension of the work begun by Leggette. It is demonstrated in Part I that by 
multiplying the cumulative departures from progressive averages by a predetermined coefficient 
and adding each in turn to a constant, there may be determined values of the ‘ “^effective average 
rate of precipitation.” In the absence of large variations from year to year in the ratios to pre- 
cipitation of the surface -rxmoff and of the loss due to evaporation and transpiration, a graph of the 
“effective average rate of precipitation” theoretically should conform with a plotting of annual 
mean water-levels. By this means an arbitrary step in scale adjustment is eliminated. Further, 
it is shown that the length of the optimum period of adjustment (in this case empirically found to 
be 25 years) is related to certain physical characteristics of the flow- system, most of which are 
quite readily determined. 

Long Island, the scene of many early groimd-water studies in the United States, affords an ex- 
cellent field for investigation. That part of the Island that lies in Nassau County and western Suf- 
folk County may be considered to represent part of an infinite strip of land of uniform width. The 
parallel boundaries of this infinite -strip aquifer are maintained virtually at constant zero-head, 
that is, at sea-level. In developing the theory it will be convenient to make further simplifying 
assumptions in regard to the h 3 qDOthetical aquifer substituted for the actual island. Discussion of 
the practical significance of these assumptions is reserved for Part n, where the effect of the 
physiography and geology on water-levels is considered. 

Recharge from precipitation is assurhed uniform over the width of the strip. Assuming now 
that the effective thickness of saturated beds below sea-level is great compared to the maximxim 
height of the water-table above sea-level so that the total thickness may be regarded as constant, 
the equation for the steady-state profile of the water-table may be expressed in the simplest 
terms. The problem is a two-dimensional one. Take as origin of coordinates a point at sea-level 
on one edge of the strip (see Fig. 1-A). Let the width of the strip be 2a. Then x = 0 and x = 2a 
represent the boundaries of the strip. Let hQ be the ordinate of the water-table profile, or the 
head, which is a function of x. Because of symmetry of the profile, at x = a there is no lateral 
flow and the slope of the water-table is zero. At any point x the rate of lateral flow per unit 
normal thickness is 


q = -Km(dhQ/dx) = W(x - a) (1) 

where m is the constant thickness of saturated beds, K is the transmission-constant or “coefficient 
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of permeability” (discharge per imit-area per xmit hydraulic gradient), and W is t^e rate of accre- 
tion to the water-table of water derived from precipitation. W has the dimensions of a velocity. 

Integrating equation (1), after multiplying through by dx 

T-hg = W(ax - x^/2) + C 

Here T, called “transmissibility,” is the product of the transmission-constant, K, and the thick- 
ness, m. It is seen that to satisfy the condition h^ = 0 at x = 0, the constant of integration, C, 
must be zero. Therefore, the equation of the steady- state profile is 

ho = (W/T) (ax - yfi/1) = (a%/2T) (2x/a - x^/a^) (2) 





Fig. 1--States of flow in hypothetical aquifer of large thickness with 
uniform accretion from precipitation: (a) Steady-state profile; 

(b) decay of profile in non- steady state 


This is the particular solution of the fundamental differential equation 

(d^ho/dx^) = -(W/T) 

obtained by differentiating equation (1), with boundary - c onditions ho = 0 at x 
(dhp/dx) = 0 at X = a. 

2 2 

The parabola (2x/a - x /a ) may be developed in a sine-series, whence 
ho = (a^W/2T) (32 /tt^) [sin(7Tx/2a) + (l/27)sin(3 7Tx/2a) 

+ (l/125)sin(5 7T x/2a) + . . ,] 

Consider now the decay of this steady-state profile in the absence of recharge, W instantane- 
ously becoming zero. Let h be the ordinate of the profile in the non-steady state, and let S be the 
“coefficient of storage,” in this case equal to the specific yield. The differential equation [5] to 
be satisfied is 

(3h/at) = (T/S) O^h/ax^) (5) 

with conditions h = 0 at x = 0 or at x = 2a for all values of t and h = ho(x) at t = 0. 

The particular solution sought here is 

h = (a^W/2T) (32 / tt^) [exp(- 7r^Tt/4a^S)-sin(7r x/2a) + 

+ (l/27)exp(-9 7T^Tt/4a^S)-sin(3 7T x/2a) + (6) 

+ (l/125)exp(-25 7T^Tt/4a^S)-sin(5 7 T x/2a) + . . .] 

Instead of following the decay of the water-table profile from h(x) = ho(x) to h(x) = 0, it will 
be desired to follow the decay from the initial steady state, corresponding to a constant rate of 


( 3 ) 

= 0 or at X = 2a, and 


(4) 
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accretion, W, to a second steady state, corresponding to a rate of accretion differing from the 
first by a small amount AW. Then the non-steady profile of the water-table is given by 

h' = (1 + AW/W)ho - (AW/W)h = Lq + (AW/W) (h^ - h) (7) 

where ho is given by equation (4) or equation (2), and where h is given by equation (6). Equation 
(7) is the particular solution of the equation 

{dh'/dx^) = -(W + AW)/T + (S/T) (3h'/at) (8) 

with boundary-conditions h' = 0 at x = 0 or at x = 2a for all values of t, and h' = ho for t = 0, satis- 
fying also the condition h' = (1 + AW/W)ho for t ->-oo. 

Because equation (8) is linear, it is permissible to superpose non-steady states described by 
equation (7). Thus theoretically it is possible in the limit to analyze the effect of rates of accre- 
tion which vary continuously. However, in correlating actual precipitation- records with water- 
level data, it will be desirable to consider the rate of accretion as varying discontinuously. Also, 
it will be foxmd permissible in dealing with Long Island to assume as a first approximation that the 
rate of accretion is constant for each of the successive (yearly) periods and proportional to the 
average rate of precipitation during that period. In other words, Wj^~Rjj, where Rj^ is the average 
rate of precipitation during the n^^ period (n^^ year). 

To simplify further the analysis which follows, neglect all terms after the first in the series 
of equation (6). Then 

h = h^-exp(- 7T ^Tt/4a^S) (9) 


and 


ho = (a^W/2T) (32/ 7T^) sin(7Tx/2a) (10) 

2 2 

In equation (9) put (4a S/ tt T) = tg and expand exp(-t/to) in a series 

(h/h„) = 1 - (t/to) + (tVato^) - (tVeto®) + (tV24t„^) - . . (11) 

It Will be foxmd expedient to replace this exponential curve by the parabola 

(h/ho) = 1 - (t/to) + (t^4to^) (12) 

terminating at t = 2to (see Fig. 1-B). 

The non- steady flow of water in the ideal system postulated above has been analyzed under 
assumptions other than those taken- -assumptions which fit the particular case more closely per- 
haps. For instance, it may be considered (1) that the sole of the aquifer is not impermeable as 
tacitly assumed above or (2) that m is not large compared to the average ho, that the total thick- 
ness varies appreciably with x [6]. Such assumptions lead to equations more complex than, but 
sufficiently near in form to, the foregoing relations to allow, for the present analysis, that they be 
approximated by equation (12). It is thought that its greater instructional value justifies treatment 
of the simplest case to the exclusion of the others at this point. Moreover, regardless of the de- 
tails of the theory there is inevitably an empirical adjustment to be made to fit it to the factual 
data. 

Rewriting equation (12) to express h as a function of t 

h(t) = h„[i . (t/y + (t2/4t„^)] (13) 

This curve may be employed to determine values of the “effective average rate of precipita- 
tion“ from past records of precipitation. By the “effective average rate of precipitation^’ at any 
given time is meant that rate of precipitation which, had it bqen maintained uninterruptedly through- 
out the past, would have been effective in producing the same water-table profile as actually existed 
at that particular time. It is assumed, as before stated, that over a sxifficiently long period the rate 
of accretion to the water-table is proportional to the rate of precipitation (W-~-R). 

Designating the effective average rate of precipitation at time t by R 
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R = R(t) dRt_ (14) 

where <|> is given by 


^ = 1 - [h(r)/hj 

(t - r) being the time of occurrence of precipitation- rate R (see Fig. 2). 

By equation (13) 

= (r/to)- (r^/4tQ^) (15) 

Hence 

— P ^^0 2 2 

R = R(t-2t„) [(-r/y - /«o )] oR/3^)(t.^dT) (16) 

Ordinarily the limits of integration here would be zero and infinity, that is, if the function 
were derived from the exponential relation of equation (9). However, for practical considerations 
the period of integration is purposely terminated, and there appears in equation (16) the constant 

Usually R is not a known continuous function and consequently (8R/9r) is not expressible in 
integrable form. Therefore, it is necessary to proceed by analyzing finite differences. The inter- 
val 2to is divided into k periods of equal length At (see Fig. 3), r is set equal to iAt, and t = nAt. 
Then dR^^_^) = -(dR/0r) is replaced by 

^R(n-i)= Vi+l)- Vi) 

is replaced by a series of coefficients having values less tlian unity as in equation (18). 

= [2ik - i(i - l)]/[k(k + 1)] (18) 



Fig. 2 --Integration of slope of precipitation-graph with $ -function to 
determine effective average rate of precipitation at time t 
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To confirm the equivalence of 0(^.1) and make the substitutions k = (2t^/At) and i = (r/At) 
and pass to the limit u > \ / 

Lim (as At -j. 0)C(^_.^ = (2 r 2t^ - rh/{2tf = (r/y - {r^/Ai^) = $ 

Replacing the integral sign of equation (16) by a summation sign with properly chosen limits 

(k-1) 

^ = %-k+l) + (^(n-i+1) - ^(n-i)) (19) 

Observing that C(n_k) = 1 and that = 0, it is seen that 

_ k 

^ " i?i “ ^(n-i+1)) (20) 



This is equivalent to the integral relation R = fR- . 

Equation (20) may be simplified further. For, as is easily demonstrated, 


Then 

^(n-i) " ^(n-i+1) - 2(k - i + l)/k(k + 1) 
_ k 

(21) 


^ " i-1 ^(""-1+1) ■ + i)A(k + 1) 

(22) 


This provides a convenient method of determining the effective average rate of precipitation 
for small values of k, though the computations become rather laborious for large values of k. 
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However, there is a still shorter method than that given by equation (22), For, as will be shown, 
the effective average rate of precipitation may be determined also simply by cumulating the de- 
partures from progressive averages of the rate of precipitation, multiplying the values thus deter- 
mined by a known rational coefficient dependent on k, and adding the appropriate initial value of ef- 
fective average rate of precipitation. 


The average rate of precipitation for k years preceding the m^^ year is (1/k) i; R/ •>,. The 

i=l ^ 
th 

departure of the rate of precipitation during the m year from the average for the k preceding 
k 

years is R^^ - (1/k) R(m-j)' 


Cumulating such departures from progressive averages, beginning with the p^ year, one gets 


(1/k) 2 [k ■ - r )] 

m=p+l j=l ^ 

It is now required to show that 

^ = [2/(k + 1)] (1/k) Z [k • R^ - Z (23) 

m=p + l J=1 ' 

where [2/(k + 1)] is the “rational coefficient’^ referred to previously. Rp is “the appropriate ini- 
tial value of effective average rate of precipitation,” which may be computed from equation (22) by 
substituting p for n 

k 

Rp = Z R(p_i+i) ' 2(k - i + l)/k(k + 1) (24) 

i=l 

After transposing the Rp to the left-hand side of equation (23), substitute for R^^ the expression 
given by equation (22), and for Rp the expression given by equation (24). Then cancel out the 
[2/k(k + 1)]. This leaves to be demonstrated that 

k n k 

^^(n-i+1) “ ^(p-i+1)^ (k - i + 1) = Z [k- R^^ - Z R(m_j)] 

1=1 > / m=p+l j=l ^ 

The right-hand member is equal to 


S t (Rm- R(m-i)) = £ Z (R„ - 

m=p+l]=l 1=1 m=p+l 

k 

which may be expanded to give Z (R, .x - R. . + R, - R, 

i=l (P+1) (p-i+1) (p+2) (p-i+2) 

• • • + R(n-l) - ^(n-i-1) + ^n ‘ ^(n-i)^' 

Assigning successive values to i 

(Rn - Rp) 4 - 

(R^ - Rp + R(n-l) - R(p-l)) + 

(Rn - Rp + R(n_i) - R(p_i) + R(n_2) “ %-2)^ + . . . + 

(Rn - Rp + R(n_i) - ^(p-l) + • ■ ■ + R(n-k+2) ‘ ^(p-k+2)^ + 

+ *^(n-l) - R{p-1) + • ■ ■ + R(n-k+2) ‘ R(p-k+2) + R(n-k+l) ‘ ®(p-k+l> 


Combining like terms 
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Table l-- Ten-year effective average rate of precipitation at the end 
of the year 1909, computed first by means of equation (19) and 
then by means of equation (22), based on records of 
precipitation at the Battery in New York City 


Year 

(n-i) 

(i) 

Rn-i 

Long method 
[using eq. (19)] 

[^n-i+1 " Rn-i)^n-i 

Short method 
[using eq. (22)1 
(Rn-i+l)(ll-i)/55 

1899 

10(=k) 

in/yr 

5.28(54/55) = 5.18 

41.78(1/55' 

) = 0.76 

1900 

9 

41.78 

47.06 2/55 

) = 1.71 

1901 

8 

47.06 

0.01 52/55) = 0.01 

47.07(3/55’ 

) = 2.57 

1902 

7 

47.07 

1.53(49/55 = 1.36 

48.60 4/55' 

) = 3.53 

1903 

6 

48.60 

-7.03(45/55 = -5.75 

41.57(5/55' 

) = 3.78 

1904 

5 

41.57 

2.91(40/55) = 2.12 

44.48(6/55' 

) = 4.85 

1905 

4 

44.48 

-2.66 34/55) = -1.64 

41.82(7/55' 

) = 5.32 

1906 

3 

41.82 

3.46(27/55 = 1.70 

45.28(8/55' 

) = 6.59 

1907 

2 

45.28 

-3.85(19/55 = -1.33 

41.43(9/55* 

) = 6.78 

1908 

1909 

1 

0 

41.43 

41.55 

0.12(10/55) = 0.02 

41.55(10/55)= 7.56 


2 ...... +1.67 43.45 

Rq9 = +1.67 + 41.78 = 43.45 = Rq9 


Table 2- - Ten-year effective average rate of precipitation at the end of the 
year 1900 computed by means of equation (24), and at the end of 
each succeeding year up to 1909 computed by means of 
equation (231. based on records of precipitation 
at the Battery in New York City 


Year 

(p-i) 

(i) 

^p-i 

Using eq. (24) 
(Rp.i^l)(ll-i)/55 




in/yr 







1890 

10(=k) 





41.441 

(1/55' 

) = 

0.75 

1891 

9 


41.44 



38.901 

2/55' 

) = 

1.42 

1892 

8 


38.90 



53.01(3/55' 

) = 

2.89 

1893 

7 


53.01 



44.17(4/55' 

) = 

3.21 

1894 

6 


44.17 



35.73(5/55* 

) = 

3.25 

1895 

5 


35.73 



37.99(6/55! 

I = 

4.15 

1896 

4 


37.99 



44.27(7/55' 

) = 

5.63 

1897 

3 


44.27 



45.121 

(8/55 

) = 

6.56 

1898 

2 


45.12 



42.061 

9/55* 

) = 

6.88 

1899 

1 


42.06 



41.78( 

:i0/55)=. 

7.60 

2 









42.34 = Rqq 

Year 

(n) 

(i) 

Rate of 
precipi- 
tation 

10-yr 

progres- 

sive 

average 

Depar- 

ture 

from 

prog. 

average 

Cumu- 

lative 

depar- 

ture 

Cum. 

dep. 

(2A1) 

Effective 
average 
rate of 
precipi;::_ 
tation, (Rj^) 


1900 

0 

in/yr 

41.78 





in/yr 

42.34 

1901 


47.06 

42.45 

4.61 

4.61 

0.84 

43.18 

1902 


47.07 

43.01 

4.06 

8.67 

1.58 

43.92 

1903 


48.60 

43.83 

4.77 

13.44 

2.44 

44.78 

1904 


41.57 

43.38 

-1.81 

11.63 

2.11 

44.45 

1905 


44.48 

43.13 

1.35 

12.98 

2.36 

44.70 

1906 


41.82 

44.00 

-2.18 

10.80 

1.96 

44.30 

1907 


45.28 

44.38 

0.90 

11.70 

2.13 

44.47 

1908 


41.43 

44.48 

-3.05 

8.65 

1.57 

43.91 

1909 


41.55 

44.11 

-2.56 

6.09 

1.11 

43.45 
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k(R„ - Rp) + (k - 1) (R(n_i) - R(p.i)) + (k - 2) (R(^_2) - R^p_2P + . . . 

+ ^(^(n-k+2) " ^(p-k+2)^ ^^(n-k+1) " %-k+l)^ 

But this is identical to 

k 

(k - i + 1) - R(p_i4.i)) 

thus proving the proposition at hand. 

The different procedures that may be followed in computing values of the effective average 
rate of precipitation are illustrated in Tables 1 and 2 and in Figure 3, using data obtained by the 
United States Weather Bureau at the Battery in New York City. A value of ten is assigned to k. 
Table 1 gives two equivalent methods of computing the ten-year effective average rate of precipita- 
tion at the end of the year 1909. First the ‘Tong method'' is given, following equation (19), in which 
finite increments in the rate of precipitation are multiplied by the appropriate coefficients 
Cn-i = [lOi - i(i - l)/2]/55 and then summed up and added to the value of R for the (n - 9)^“ year 
(the year 1900). The “short method," following equation (22), involves merely summing up 
products of rates of precipitation and appropriate differences between successive coefficients: 
(C(„_i) -C(„.ul))=(ll-l)/55, 

In Table 2 the ten-year effective average rate of precipitation at the end of the year 1900 is 
computed by means of equation (24), “short method". Values of the ten-year effective average 
rate of precipitation at the end of succeeding years are computed by means of equation (23). First, 
values of ten-year progressive average precipitation are determined. Then, the departures of the 
actual annual rates of precipitation from the progressive averages are determined and cumulated. 
The cumulated departures are then multiplied by [2/(k + 1)] = (2/11) and added to the effective 
average for 1900 (Rqq). 

It is seen that the effective average for 1909, determined by three different procedures, is 
43.45 inches per year. Computations shown in Tables 1 and 2 were carried out on a calculating 
machine to one place beyond what is tabulated. Accordingly, minor adjustments were made in 
rounding off the figures in the last columns of both Tables. 
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Part II- -Correlation of data 

Abstract - -Records of rainfall at the Battery in New York City and at Setauket, Long Island, 
are analyzed by the procedure outlined in Part I. In determining progressive averages, periods of 
10, 15, 20, 25, 30, and 40 years are used. Values of the effective average rate of precipitation are 
correlated with composite average water-levels for the periods 1912-18 and 1932-41 based on 
records of 14 wells in Nassau County and western Siiffolk Coimty. The best correlation is ob- 
tained with the 25 -year effective average, a fair agreement being shown with annual mean water- 
levels. Analysis and correlation of the same data on a 300-month basis permit an approximate 
evaluation of the relative magnitude of the loss due to transpiration and evaporation. Knowledge 
of the approximate value of this fundamental time-constant (25 years or 300 months) permits a 
crude checking of the probable values of other factors such as the permeability and specific yield 
of the beds and their effective thickness, and the average annual rate of accretion to the water- 
table. 

For the purpose of comparison, correlations of water-levels are also made using (1) pro- 
gressive averages of precipitation and (2) cumulative departures from fixed averages- -two methods 
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of presentation now commonly used. 

The effective average,'' determined by cumulating departures from progressive averages, 
is suggested as a tool for the analysis also of long-term records of “ground-water flow" in 
streams. 

U. S. Geological Survey, 

Jamaica, New York 



CORRELATION OF GROUND-WATER LEVELS AND PRECIPITATION 
ON LONG ISLAND, NEW YORK 

C. E. Jacob 

(Published with the approval of the Director, U. S. Geological Survey) 

PART n- -CORRELATION OF DATA 

A brief though concise statement of the history of ground-water studies on Long Island, begin- 
ning with the early water-level observations in Brooklyn by STODDARD in 1854, was given by 
THOMPSON [see 7 of ‘'Referenced" at end of paper]. These and other early data were considered 
later by LEGGETTE [8]. He evaluated them by means of a graph of the cumulative departure of 
precipitation. More recent studies by LEGGETTE [9] and by the writer have lead to the procedure 
outlined in Part I of this paper [10], which was founded upon an empirical approach suggested by 
LEGGETTE and was later justified by analysis based on the theory of BOUSSINESQ [11]. 

Figure 1 is a map of the central part of Long Island showing contours on the water-table as of 
May, 1943, and also the locations of 14 New York City test-wells upon whose records the present 
study is based. The contours are based upon data from many more wells than are actually shown 
on the map. The 14 wells whose water-levels are analyzed in this study were first measured by 
the Department of Water Supply, Gas, and Electricity of the City of New York. They were chosen 
because of the length and continuity of their records and because together they adequately cover 
the range of water-table elevations in the area with more or less uniform distribution. Other wells 
with records of eqxml or greater length were eliminated because of less favorable situation with 
respect to centers of pumpage, points of natural discharge, etc., or because they duplicated records 
of one of the 14 wells. 



configuration of main ground-water table in May, 1943, and location of shallow test-wells 

9 
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As shown by the last two columns of 
Table 1, the approximate depth to water 
in the wells is on the average about half 
the elevation of the water -surface in 
the same wells though there is con- 
siderable variation from well to well. 
Although the depth of the water-table 
below the ground -surface may be a de- 
termining factor in the relation between 
rainfall and changes in ground -water 
level during and following individual 
storms, it has much less bearing upon 
long-term or secular variations of 
water-level such as are considered in 
this study. 

The principal centers of net with- 
drawal from the shallow water -producing 
beds in the area under consideration are 
the infiltration-galleries and shallow- 
well stations of the City of New York, in 
the southern part of Nassau County. 
During the years 1912 to 1918 the aver- 
age rate of withdrawal from this system 
was about 38 million gallons per day 
(mgd), whereas from 1932 to 1941 the 
average rate of withdrawal was only 
about 23 mgd. However, the observation- 
wells on whose records the present study 
is based are all sxifficiently distant from 
these centers of pumpage that the effect 
on their water-levels may be neglected. 

Table 2 gives “composite average 
water-levels'’ and “effective average 
rates of precipitation” by years for the 
periods 1912 to 1918 and 1932 to 1941. 

The yearly values of composite average 
water-level were computed from annual 
mean water-levels in 14 or fewer wells. 
The number of wells for which there is 
s uff icient record to determine an annual 
mean water-level for a given year is 
found in the third column of Table 2. For 
the years 1914 and 1937 to 1942, mea- 
surements were frequent enough and 
were so distributed that annual mean 
water-levels could be obtained for all 14 
wells. For anyone of those years the 
composite average water-level is merely 
the average of the 14 values of mean 
annual water-level. The composite aver- 
age water-level for anyone of the other 
years was determined by comparing the 
average water-level in the wells for 
which records are available for that 
given year with the average water-level 
in the same wells during 1914 and 1937 
to 1942. That is to say, if the number of 
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Table 2- - Comparison between composite average water-level ajid effective 
average rate of precipitation on Long Island. New York 


Year 

(1) 

Elevation 
composite 
av. water- 
level, h 

(2) 

No. 

of 

wells 

( 3 ) 

Effective average rate of preciflitation 
Battery, New York 

25- _ 
year, R 

( 4 ) 

( 5 ) 




ft 


in/yr 

in/yr 

in/yr 

in/yr 

1912 

48.53 

2 

42.53 

-0.51 

42.66 

-0.43 

1913 

47.65 

11 

42.55 

+ 0.39 

42.43 

+ 0.19 

1914 

47.77 

14 

41.76 

-0.62 

42.23 

-0.16 

1915 

47.74 

7 

41.62 

-0.75 

41.79 

-0.62 

1916 

46.96 

4 

40.92 

-0.76 

41.32 

-0.37 

1917 

45.85 

3 

40.72 

+ 0.12 

40.73 

+ 0.07 

1918 

45.10 

2 

40.07 

+ 0.14 

40.37 

+ 0.42 

1932 

43.34 

6 

40.01 

+ 1.83 

40.02 

+ 1.78 

1933 

44.45 

9 

40.76 

+ 1.57 

40.41 

+ 1.11 

1934 

45.86 

6 

41.13 

+ 0.58 

40.99 

+ 0.36 

1935 

46.41 

10 

40.53 

-0.65 

40.86 

-0.34 

1936 

46.13 

3 

40.99 

+ 0.15 

40.78 

-0.15 

1937 

46.84 

14 

41.58 

+ 0.10 

41.34 

-0.23 

1938 

47.54 

14 

42.01 

-0.11 

41.78 

-0.43 

1939 

48.66 

14 

41.60 

-1.69 

42.07 

-1.23 

1940 

47.46 

14 

41.91 

-0.14 

41.84 

-0.29 

1941 

46.76 

14 

41.52 

+ 0.12 

41.73 

+ 0.29 

Mean 

46.65 

... 

41.31 

0.602 

41.37 

0.498 


= [(46.65 ft) R/(41.31 in/yr) - h]. = [(46.65 ft) R'/(41.37 in/yr) - E]. 


wells for which records are available is n, the composite average water-level for a given year is 
determined by 


(avg water-level in all 14 wells for 1914, 1937-42) 
(avg water-level in n wells for 1914, 1937-42) 


X (average water-level 
in n wells for the 
given year) 


Values of composite average water-level were determined also month by month. These val- 
ues are plotted twice in Figure 2 --in the upper blocks of that diagram to the same scale as the 
water-level data for the several wells, and in the lower blocks of the diagram to five times that 
scale. The composite average water-level is seen to follow closely the trend of water-level in 
wells N-1263 and N-1259. It is believed that this “composed’' or “composite” average water-level 
graph would practically duplicate the mean water-level graph of the 14 wells that would have been 
obtained if all those wells had been measured each month during the period of record. The method 
by which this continuous though artificial record was composed was as follows; 

The monthly composite average water-levels for 1942 are simply arithmetic means. For the 
years 1936 to 1941 the monthly composite average water-levels were determined by multiplying 
the corresponding average water-levels for eight of the 14 wells by the constant 1.044. As indicated 
in Table 1, this constant was determined by taking the ratio of the average for all 14 wells of the 
mean water-levels for 1942 to the corresponding average for the eight wells.' These eight wells are 
those that have continuous records of monthly readings during the period 1936-41. For the years 
before 1936 the monthly composite average water-levels were determined by means of the coeffi- 
cients (C) given in the third column of Table 1, the derivation of which is explained in a footnote to 
that Table. For example, if at the end of a given month water-level measurements were made in 
only four wells, the composite average water-level for that month was determined by multiplying 
the elevation of the water-level in each of the four wells by the appropriate coefficient, summing 
such products, and dividing the sum by four. 
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Values of the 25-year effective average rate of precipitation plotted in Figure 3 were computed 
by the method outlined in Table 2 of Part I of this paper [12] from data obtained by the United States 
Weather Bureau at the Battery in New York, New York. Values for the years 1912 to 1918 and 1932 
to 1941 are tabulated in Table 2 (Part H). Also there are tabulated in Column (6) of that Table year- 
ly mean values of the 3 00 -month effective average rate of precipitation for the same years computed 
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Fig. 4 — Graph of average relative error of linear correlation 
between effective average rates of precipitation at the Battery 
and composite average water-level on Central Long Island 


by the same method. Values of 300-month effective average rate of precipitation, from which 
these yearly means were determined, are plotted in the lower blocks of Figure 2, not only for the 
Battery but also for the station at Setauket, Long Island. Also there is given the average for the 
two stations. 

Effective average rates of precipitation were also determined on the basis of 10-, 15-, 20-, 

30-, and 40-year periods of adjustment from the same rainfall-data. It was foimd by trial that the 
best correlation with the composite average water-level for both periods of record on a yearly basis 
was obtained with the 25 -year effective average. This may be seen from Figure 4, in which are 
plotted values of average relative error of the correlation against the period of adjustment, kAt, 

The meaning of “average relative error” is illustrated in Table 2. The average elevation of the 
composite average water-level (h) for the 17 years of record is 46.65 feet. The average value of 
the 25-year effective average rate of precipitation (R) for the same 17 years is 41.31 inches per 
year. If the basic assumptions o^the theory in Part I were met, there should be continuously a 
direct proportion between h_and R of constant value: (46.65 ft)/(41.31 in/yr). On this basis one 
may multiply the values of R by this ratio and subtract the corresponding values of H to obtain 
the residuals A. The average of the absolute values of the residuals is in this case 0.602 foot. 
Referring this to the average value of h, the relative value of the average residual, or the “average 
relative error”, is found to be (0.602 ft)/(46.65 ft) = 0.0129. 

Actually, correlations were first effected with the water-level records of individual wells of 
the group of 14. Then, to simplify the work, and in a sense to summarize those correlations, the 
composite average water-level was adopted. It was found, as was to be expected, that the corre- 
lation could be improved further by using instead of the values of 25 -year effective average given 
in Column (4) of Table 2 the yearly means of the 300-month effective average, given in Column (6) 
of that Table. It should be pointed out that the 25 -month effective average and the 300 -month effect- 
ive average involve the same period of adjustment (25 years = 300 months) and that they are based 
on the same basic data, but that the latter divides the period into 12 times as many intervals, A 
comparison between the two may be effected by connecting year-end values of the 300-month effect- 
ive averages plotted in Figure 2 with straight-line segments, as the year -end values of necessity 
coincide with the corresponding values of the 25 -year effective average. 
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The yearly mean values of the 300-month effective average are indicated in Figure 2 by 
bars extending over each year's record. The average of these for the 17 years of record is shown 
in Table 2 as 41,37 feet. Multiplying successive values of R' by the ratio (46.65 ft)/(41.37 in/yr) = 
1.128 ft/(in/yr) and subtracting corresponding values of 1, one gets the residuals A', which average 
0.498 foot in absolute value. Here, then, the average relative error is 0.0107. This value is plotted 
on the 25 -year line in Figure 4 as a bold open circle in contrast to the light open circle on the same 
line at 1.29 per cent representing the average relative error of the correlation between the com- 
posite average water-level and the 25 -year effective average rate of precipitation. 

The smooth curve drawn through the points for the correlations based on the 17 years of record 
(1912 to 1918 and 1932 to 1941) shows a minimum at kAt = 26 years. Similarly the smooth curves 
for correlations based on the period 1912 to 1918 and on the period 1932 to 1941 both show minima 
at about 24 or 25 years. 

There is plotted in Figure 3, as a light line on the upper block of that diagram, a curve of the 
25 -year progressive average rate of precipitation. A given point on that curve represents the aver- 
age rate of precipitation at the Battery for the preceding 25 years. That curve was also correlated 
with the composite average water-level over the same 17 years of record in the manner described 
above for the effective average precipitation-rates. The average relative error was found to be 
0.0264, as plotted in Figure 4. 

Correlations were also effected between the composite average water-level and effective 
average rates of precipitation at Setauket, Long Island. Inasmuch as the record at this station only 
goes back to 1886, or only 26 years beyond 1912, it is not possible to compute 30- and 40-year 
effective averages for the entire period of water-level records. Data for the 10-, 15-, 20-, and 
25-year effective averages suggest that the optimum period of adjustment is again about 25 years, 
with an average relative error of about 0,016. 

An effort was made empirically to improve upon the linear correlation between the composite 
average water-level and the yearly means of the 300 -month effective average pr e cipitation - r ate by 
adding an arbitrary constant to values of the former before determining their ratio to values of the 
latter. The best correlation was obtained by adding 15 feet to the composite average water-level. 
The ratio then became (61.65 ft)/(41.37 in/yr) = 1.490 ft/(in/yr). Values of composite average 
water-level for 1912 to 1918 and 1932 to 1942 adjusted in this manner are plotted in Figure 3 as 
open circles falling near the composite average water-level curve on the upper block of that dia- 
gram. Monthly values of effective average precipitation-rate adjusted in the converse manner are 
plotted on the lower block of Figure 2. The resulting curve runs approximately through the mean 
position of the curve of the composite average water-level. 

The fact that the best correlation is obtained not by a simple proportion but that the addition of 
a constant is required may be attributed to the failure of conditions on Long Island to satisfy com-^ 
pletely two of the basic assumptions of the theory. The first, and perhaps more important, is that 
“the effective thickness of saturated beds below sea-level is great compared to the maximum 
height of the water-table above sea-level so that the total thickness may be regarded as constant" 
[Part I, p. 565]. The second is that “the rate of accretion is constant for each of the successive 
periods and proportional to the average rate of precipitation during that period" [Parti, p. 567]. 

It was pointed out in Part I [p. 567] that a modified theory taking account of the lateral variation of 
thickness of the saturated beds leads to equations sufficiently near in form to those of the simple 
theory that the approximations of the latter may be accepted, especially in view of the empirical 
adjustment that is required. Taking another view in justification of the procedure followed, the 
best linear correlation between the compdsite average water-level and the effective average rate 
of precipitation was sought in order to improve upon extrapolations of water-level records into the 
past and yet to insure reasonable overall conformity between current water-level and rainfall -data 
as a basis upon which to work in analyzing short-term water-level fluctuations. Having accounted 
for the secular variations of ground-water level, we should be better prepared to proceed with the 
analysis, for example, of the annual phreatic cycle. 

Referring again to Figure 2, and comparing the curve of composite average water-level with 
the curve for the 300-month effective average precipitation-rate adjusted linearly to the same 
scale, it may be remarked that although over a period of years the two curves agree on the average, 
in detail they do not agree. But the assumptions upon which the theory is based, especially the one 
regarding the proportionality between rates of precipitation and rates of accretion, are such that a 
closer agreement than that existing between these two curves is not to be expected. It was stated 
in Part I [p. 565] that “in the absence of large variations from year to year in the ratios to 
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precipitation of the surface -runoff and of the loss due to evaporation and transpiration, a graph of 
the effective average rate of precipitation theoretically should conform with a plotting of annual 
mean water-levels." Lack of conformity between the two curves on an annual basis, as plotted in 
Figure 3, may be attributed, therefore, to variation in the fractions of the yearly precipitation that 
are disposed of as runoff and as evaporation and transpiration. Similarly, the more evident lack of 
conformity between the same two curves plotted on a monthly basis in Figure 2 may be attributed 
to the greater month -to -month and seasonal variations in the same factors. Major causes of dis- 
crepancy are the occurrence of successive storms of high intensity, as illustrated by the record 
for 1938 and 1939 [13], and the prolongation of periods of drought, as for example in the late sum- 
mer and fall of 1941. 

An idea of the relative magnitude of the seasonal variation in evaporation and transpiration as 
reflected in the annual cycle of the water-level may be gained by comparing the amplitudes of the 
two curves referred to, particularly for the more or less "normal" years 1912 to 1915 and 1935 to 
1937. According to the theory outlined in Part I, if, starting on a given day, for example at the 
phreatic high in the spring, all further accretion to the water-table were to cease, the water-level 
would begin to decline at a rate such that if this initial rate of decline were continued the water 
table would be reduced to sea level in the time tp = (kAt/2), which in this case is 12 or 13 years. 
Actually the rate of decline of the water-table during the late spring and early siunmer of the years 
cited was on the average about twice this theoretical maximum. Of course allowance should be made 
for the possibility that a more exact theory taking account of the non-linearity introduced by the 
variation in depth of flow would allow a higher rate of initial decline during a period of zero -accre- 
tion. Nevertheless, even with this in view it is suggested quite clearly that over a fair proportion 
of the annual cycle on Long Island the movement of moisture from the ground to the water-table is 
reversed. The rate of accretion thus becomes temporarily negative. The motivating processes 
are, of course, evaporation and transpiration. But the exact details of the mechanism of transfer 
remain to be elucidated. 

In Figure 3 there is plotted a curve giving the cumulative departure of precipitation from the 
average rate at the Battery for the period 1826 to 1942 [14] . This may be compared with the graphs 
of the 25 -year progressive average pr e cipitation - r ate and the 25-year effective average precipita- 
tion-rate. Of the three curves, that of the 25 -year effective average correlates best with the com- 
posite average water-levels plotted as open circles. Unfortunately the water-level measurements 
only cover a range of secular fluctuation of about four feet in 47 feet, or less than ten per cent. 
However, we do have reports by the natives on the level of Lake Ronkonkoma going back to 1857. 

The statements given in Figure 3 are taken from the report of the BURR-HERENG-FREEMAN 
Commission [15] referred to by LEGGETTE [16]. Making due allowance for the fallibility of mem- 
ory, those reports check fairly well the probable trend of water-levels in those years inferred from 
the 25-year effective average precipitation-rate. The degree to which the level of Lake Ronkonkoma 
follows the trend of ground- water levels in the area to the west is seen in Figure 2, where a month - 
by-month record of the lake-level is given for recent years. 

Having determined empirically the value of 1 q, we are enabled crudely to check certain other 
factors to which it is related. In Part I, following equation (10), tg is designated as an abbreviation 
for (4a2s/7T^T), where a is the half -width of the hypQthetically iniinite strip, T is the transmissi- 
bility of the aquifer, which is assumed to be of great thickness, and S is the coefficient of storage, 
in this case an "effective porosity". From equation (2) of Part I, when x = a, ho = (a2w/2T). This 
gives the elevation of the water-table at the ground-water divide in terms of a, T, and W, the latter 
being the effective average rate of accretion to the water-table. Substituting in the foregoing ex- 
pression for to the value of (a2/T) given by the last equation above 

to = (4S/7r2)(2ho/W), or W = (SShoA^to) 

Now, it was found that to - 12.5 years = 4,600 days. On Section A-A' (Fig. 1) ho = 85 feet 
(for X = a). And S may be taken as approximately 0.40. Thus W = (8 x 0.40 x 85 ft)/(7T x 4,600 
days) = 0.0060 foot/ day. 

The effective average precipitation-rate is about 43 in/yr (average for the Battery and Setauket), 
or R = 0.0098 foot/day. This would indicate that the average rate of accretion to the water-table is 
about 6p per cent of the average rate of precipitation. An earlier estimate placed this ratio for the 
island as a whole at about 44 per cent [17]. However, there is sufficient leeway in both estimates 
that the difference need not be considered significant. 
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The transmissibility of an equivalent aquifer of large thickness may now be determined: 

T = (a2w/2ho) = (44,0002ft2 x 0.0060ft/day)/(2 x 85 ft) = 68,OOOft2/day 

The greatest uncertainty, perhaps, attaches to the proper value of thickness m. If it be assumed 
that the sands of Magothy age are effectively inter-connected with the water-table, then an average 
thickness of 500 feet might be taken as reasonable; whence the transmission-constant K = (T/m) = 
140 ft/ day . It should be emphasized, however, that the above value of T is only an apparent trans- 
missibility. Moreover, the equating of the two expressions for (a^/T) to determine W is not strict- 
ly valid and gives, therefore, only a rough approximation of the true average W. 

The concept df effective average rate of precipitation might find useful application in the analy- 
sis of long-term records of “ground- water flow" in streams. Indeed, it would seem that a step 
toward bridging the gap between our understanding of infiltration and of ground-water runoff might 
be taken by setting up in the field a system of observation-wells and gaging -stations in as near an 
ideal area as possible --one that is not unduly complicated by effects of frost or snow- cover. If an 
area could be foimd in a humid climate where a nearby flat-lying bed-rock is covered more 
or less uniformly with sediments cut by parallel streams, the installation of a series of test-well 
profiles and rain-gages, with one or more pairs of gaging- stations, would provide sufficient control 
for a rational evaluation of ground-water runoff. One of the principal weaknesses of the present 
study on Long Island lies in the impracticability of measiiring directly the ground-water that is 
discharged along the northern and southern shores of the Island. In the suggested test set-up these 
parallel shores would be replaced by effluent streams. The gain in flow per unit-length of reach 
would give, then, the desired measure of ground-water discharge. 
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